Response of thin-film SQUIDs to applied fields and vortex fields: Linear SQUIDs 
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In this paper we analyze the properties of a dc SQUID when the London penetration depth A 
is larger than the superconducting film thickness d. We present equations that govern the static 
behavior for arbitrary values of A = X^/d relative to the linear dimensions of the SQUID. The 
squid's critical current 7c depends upon the effective flux <I>, the magnetic flux through a contour 
surrounding the central hole plus a term proportional to the line integral of the current density 
around this contour. While it is well known that the SQUID inductance depends upon A, we show 
here that the focusing of magnetic flux from applied flelds and vortex-generated flelds into the central 
hole of the SQUID also depends upon A. We apply this formalism to the simplest case of a linear 
SQUID of width 2w, consisting of a coplanar pair of long superconducting strips of separation 2a, 
connected by two small Josephson junctions to a superconducting current-input lead at one end and 
by a superconducting lead at the other end. The central region of this SQUID shares many properties 
with a superconducting coplanar stripline. We calculate magnetic-field and current-density profiles, 
the inductance (including both geometric and kinetic inductances), magnetic moments, and the 
effective area as a function of A/w and a/w. 

PACS numbers: 74.78.-w 
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I. INTRODUCTION 



The research in this paper has been motivated by 
several important recent developments in superconduc- 
tivity: (a) the fabrication of thin-film SQUIDs (su- 
perconducting quantum interference devices) made of 
high- Tc superconductors,^- (b) the study of noise gener- 
ated by vortices in active and passive superconducting 
devices fiiS42iiiSiSiLSi2ii& and (c) line-width reduction in 
superconducting devices to eliminate noise due to vor- 
tices trapped during cooldown in the earth's magnetic 

n.12.1.T14.1.'^.1fi 

The fact that the London penetration depth A in- 
creases as T increases and diverges at Tc is an important 
consideration for high-Tc SQUIDs operated at liquid- 
nitrogen temperature. When A is larger than the film 
thickness d, the physical length that enters the equations 
governing the spatial variation of currents and fields is 
the Pearl lengtbii A = A^/d. Accordingly, the equations 
governing the behavior of active or passive thin-film su- 
perconducting devices depend upon the ratio of A to the 
linear dimensions of the device. In particular, the equa- 
tions governing SQUIDs involve not just the magnetic 
fiux up through a contour within the SQUID, but the ef- 
fective flux <i>, which is the sum of the magnetic flux and 
a term proportional to the line integral of the current 
density around the same contour. While the effective 
flux $ is similar to London's fluxoid,^^ which is quan- 
tized in multiples of the superconducting flux quantum 
00 ~ h/2e, we show in the next section that $ is not 
quantized. We also give in Sec. II the basic equations, 
valid for any value of A, that govern the behavior of a dc 
SQUID. 



A vortex trapped in the body of the SQUID during 
cooldown through the superconducting transition tem- 
perature Tc in an ambient magnetic field generates a mag- 
netic field and a screening current that together make a 
sizable vortex-position-dependent contribution to the 
effective fiux $. If such a vortex remains fixed in po- 
sition and the temperature remains constant, this sim- 
ply produces a harmless bias in $. On the other hand, 
both vortex motion due to thermal agitation and temper- 
ature fluctuations generate corresponding fluctuations in 
<i>v and noise in the SQUID output. In Sec. II we show 
that to calculate $v, it is not necessary to calculate the 
spatial dependence of the vortex-generated fields and cur- 
rents. Instead, one may determine <I>v with the help of 
the sheet-current distribution of a circulating current in 
the absence of the vortex. 

In Sec. Ill we apply the basic equations of Sec. II 
to calculate the properties of a model linear SQUID, 
which has the basic topology of a SQUID but is greatly 
stretched along one axis, such that the central portion 
resembles a coplanar stripline. The advantage of using 
such a model is that simple analytical results can be 
derived that closely approximate the exact numerically 
calculated quantities in the appropriate limits. In ad- 
dition to calculating the field and current distributions 
for several values of A, we calculate the total inductance, 
geometric inductance, kinetic inductance, and magnetic 
moment when the SQUID carries a circulating current. 
We calculate the field and current distributions, magnetic 
moment, and effective area Acs ~ ^t/Sa when a perpen- 
dicular magnetic induction i?a is applied and the effective 
flux is focused into the SQUID. Finally, we calculate 
the field and current distributions and the magnetic mo- 
ment for the zero-fluxoid state when the junctions are 
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short-circuited and the sample remains in the state with 
$ = when a perpendicular magnetic induction i?a is 
applied. 

In Sec. IV, we present a brief summary of our results. 



II. BASIC EQUATIONS 

Our purpose in this section is to derive general equa- 
tions that govern the behavior of a dc SQUID consisting 
of thin superconducting films of thickness d less than the 
weak-field London penetration depth A, such that the 
fields and currents are governed by the two-dimensional 
screening length or Pearl length-^ A = A^/d. We calcu- 
late both the current I = I1+I2 through the SQUID [see 
Fig. 1] and the circulating currenliiS 1^ = (/2 — /i)/2 and 
describe how to calculate the critical current Ic of the 
SQUID for arbitrary values of the SQUID's inductance 
L. For this case, the contributions from line integrals of 
the current density to the effective flux in the hole cannot 
be neglected, and the kinetic inductance makes a signifi- 
cant contribution to L. When a perpendicular magnetic 
induction i3a is applied, we calculate how much magnetic 
flux is focused into the SQUID's hole; this flux also can 
be expressed in terms of the effective areaSii of the hole. 
We also show how to calculate how much magnetic flux 
generated by a vortex in the main body of the SQUID is 
focused into the hole. 

Consider a dc SQUID in the xy plane, as sketched 
in Fig. 1. We suppose that the SQUID is symmetric 
about the y axis, which lies along the centerline. The 
maximum Josephson critical current is Iq for each of the 
Josephson junctions, shown as small black squares. The 
currents up through the left and right sides of the SQUID 
can be written as Ii = 1/2 — Id and I2 = 1/2 + Id- 
When the magnitude of /, the total current through the 
SQUID, is less than the critical current 1^, the equations 
that determine I = Ii + I2 and the circulating current 
-^d = {h ~ -^i)/2 can be derived using a method similar 
to that used in Ref. |^ We begin by writing the local 
current density j in the superconductors (i.e., the main 
body of the SQUID and the counterelectrode) asi^ 



j = -(l/AioA')[A + (^o/27r)V7], 



(1) 



where A is the vector potential, 0o — h/2e is the super- 
conducting flux quantum, and 7 is the phase of the order 
parameter. The quantity inside the brackets, which is 
gauge- invariant, can be thought of as the superfluid ve- 
locity expressed in units of vector potential. From the 
point of view of the Ginzburg-Landau theory, implicit in 
the use of this London-equation approach is the assump- 
tion that the applied flelds and currents are so low that 
the magnitude of the order parameter is not significantly 
reduced from its equilibrium value in the absence of fields 
and currents. 

To obtain the SQUID equations, we integrate the vec- 
tor potential around a contour C that passes in a coun- 
terclockwise direction through both junctions, the main 




FIG. 1: / enters the main body of the SQUID from the coun- 
terelectrode below through two Josephson junctions (small 
black squares, labeled a and b) and divides into the currents 
Il — 1/2 — /d and I2 = 1/2 + Jd as shown. 



body of the SQUID, and the counterelectrode as shown in 
Fig. 1 and write the result in two ways. Since B = V x A, 
this integral yields, on the one hand, the magnetic flux 
in the z direction 



B^{x,y)dS, 



(2) 



where S is the area surrounded by the contour C and 
Bz {x, y) is the z component of the net magnetic induction 
in the plane of the SQUID produced by the sum of a 
perpendicular applied held Bg_ and the self-field Bsz{x^ y) 
generated via the Biot-Savart law by the supercurrent 
density j{x, y). On the other hand, for those portions of 
the contour lying in the superconductors, we eliminate 
the line integrals of A in favor of line integrals of j using 
Eq. (Q. We then express the line integrals of V7 in 
terms of the values of 7 at the junctions. Equating the 
two expressions for the line integral of A, we find that 
the effective flux $ in the z direction through the SQUID 
is given by 



where 



((/)o/2^)(01-02) 



<&= / B,{x,y)dS + j j 
Js Jc 



(3) 



dl, (4) 

Is Jc 

with the integration contour C now passing through 
both superconductors but excluding the junction bar- 
riers. These equations are equivalent to Eq. (8.67) in 
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Ref. 1^ The gauge-invariant phase differences across the 
junctions b and a are, respectively)^^ 

01 = 76c - 76s - (27r/0o) / A dl, (5) 

J be 
nas 

</'2 = lae - las - (27r/(/)o) / A ■ dl, (6) 

ae 

where be labels a point on the counterelectrode side of 
the junction b and bs labels the point directly across 
the insulator in the SQUID washer, and ac and as la- 
bel corresponding points for junction a. According to 
the Josephson equations?^ the junction supercurrents 
are Ii = /osin^i and I2 = /osin(/)2. In the above 
derivation we have assumed that the linear dimensions of 
the Josephson junctions are much less than the Joseph- 
son penetration depth Xjr^ and that the applied fields 
are sufficiently small that the Josephson current densi- 
ties and gauge-invariant phase differences are very nearly 
constant across the junction areas. 

The magnetic moment m = mz generated by the cur- 
rents in the SQUID is^^ 

m = i j r X jd^r. (7) 

It can be shown with the help of the London fluxoid 
quantization condition^i^ 

/ Bz{x, y)dS + ^oA^ / j dl = n0o, (8) 
Js' Jc 

where n is an integer and C is a closed contour that sur- 
rounds an area S' within the body of the SQUID, that 
if there are no vortices present (i.e., when n = 0), the 
expression for $ in Eq. Q is independent of the choice 
of contour C. Any convenient path can be chosen for 
C, provided only that the path remains in the super- 
conducting material in the body of the SQUID and the 
counterelectrode. On the other hand, when there are 
vortices in the main body of the SQUID, the quantity <& 
increases by 0o each time the contour C is moved from 
a path inside the vortex axis to one enclosing the vortex 
axis. Thus, without specifying the precise contour C, $ is 
determined only modulo 4>o. However, this is of no phys- 
ical consequence, because the gauge-invariant phases 
and (l>2 , which also enter Eq. Q , are also determined only 
modulo 27r. The final equations determining the currents 
/ and /(J are independent of the choice of contour C and 
remain valid even when vortices are present in the main 
body of the SQUID. 

When the thickness d of the SQUID is much larger 
than A, the contours C and C can be chosen to be at the 
midpoint of the thickness, where j is exponentially small, 
such that the line integrals of j can be neglected. The 
resulting equations are then the familiar ones found in 
many reference books, such as Refs. 21,24,25 26 27 28 29. 
However, we are interested here in the case for which 
d < X, such that the fields and currents are governed by 



the two-dimensional screening length or Pearl lengtbii 
A — /d. The term in Eq. ^ involving j then must 
be carefully accounted for. For this case, j is very nearly 
constant over the thickness and it is more convenient to 
deal with the sheet-current density J{x, y) — jd, such 
that Eqs. 0) and © take the formSa 

/" B^{x,y)dS + fioA [ J dl (9) 
Js Jc 

and 

/ B^{x,y)dS + noA J-dl = n4)a. (10) 
Js' Jc 

For the general case when the SQUID is subject to a 
perpendicular applied magnetic induction Sa, carries a 
current / unequally divided between the two arms, Ii = 
1/2— Id and I2 = 1/2+1^, where the circulating currenliiS 
is Id — (/2 — /i)/2, and contains a vortex at the position 
rv in the body of the SQUID, the effective flux <& in the z 
direction can be written as the sum of four independent 
contributions: 

$ = $/ + $d + 'I'f + $v, (11) 

The first term on the right-hand side of Eq. (|ll|l is that 
which would be produced by equal currents //2 in the y 
direction on the left and right sides of the SQUID shown 
in Fig. 1: 

/ Biix,y)dS + f^oA [ Ji dl, (12) 
Js Jc 

where Bj{x, y) is the z component of the self- field gener- 
ated via the Biot-Savart law by the sheet-current density 
Ji{x,y), subject to the condition that the same current 
1/2 flows through the two contacts a and b. For a sym- 
metric SQUID, Jj(x,y), the y component of Ji(x,y), is 
then an even function of x, and Ji{x, y) and Bi{x, y) are 
odd functions of x. As a result, both terms on the right- 
hand side of Eq. H12|) vanish by symmetry, and $/ = 0. 
Since V • .7/ = except at the contacts a and 6, we may 
write Ji — — (//2)V x G/, where Gj — zGj, such that 
Ji{x,y) = {I /2)zxWGi{x,y). The contours of the scalar 
stream function Gi{x,y) — const correspond to stream- 
lines of Ji{x,y), and we may choose G/ = for points 
Ti = {xi,yi) all along the inner edges of the supercon- 
ductors and G/ = 1 for points Tq — {xo,yo) all along the 
outer right edges and G/ = —1 for points To = (xo,yo) 
all along the outer left edges. 

The second term on the right-hand side of Eq. (|ll|l is 
due to the circulating currenliiS Id — {I2 — A)/2 in the 
counterclockwise direction when unequal currents flow in 
the two sides of the SQUID shown in Fig. 1: 

$d= / Bd{x,y)dS + noA ( Jd dl, (13) 
Js Jc 

where Bd{x,y) is the z component of the self-field gen- 
erated via the Biot-Savart law by the circulating sheet- 
current density Jd{x,y) when a current Id fiows through 
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contact a from the counterelectrode into the body of 
the SQUID, passes around the central hole, and flows 
through contact h back into the counterelectrode. The 
magnetic moment generated by the circulating cur- 
rent is proportional to Id, as can be seen from Eq. ((T)). 
Since V • Jd = except at the contacts a and 6, we may 
write Jd — ~-^dV x Gd, where Gd = -zGd, such that 
Jd{x,y) — IdZ X VGd(a;,y). The contours of the scalar 
stream function Gd{x,y) = const correspond to stream- 
lines of Jd(x,y), and we may choose Gd = for points 
Ti = {xi.yi) all along the inner edges of the supercon- 
ductors and Gd — 1 for points Tq — {xo, yo) a-U along the 
outer edges. Once a numerical result for $d is found, the 
result can be used to determine the inductance L of the 
SQUID via L = ^d/Id, as was done for a circular ring in 
Ref. 31. The resulting expression for L is the sum of the 
geometric and kinetic inductances. 

The third term on the right-hand side is a flux-focusing 
term due to the applied field: 

$f = / Bi{x,y)dS + fio^ [ Jfdl, (14) 
Js Jc 

where Bf{x, y) is the z component of the net magnetic in- 
duction in the plane of the SQUID produced by the sum 
of a perpendicular applied field i?a and the z component 
of the self-field Bsi{x,y) generated via the Biot-Savart 
law by the sheet-current density J{{x,y) induced in re- 
sponse to i?a, subject to the condition that no current 
flows through the junctions a and b. In other words, the 
desired fields are those that would appear in response 
to -Ba if the junctions a and b were open-circuited. Since 
V- Jf = 0, we may write Jf — —V x Gf , where Gf = zGf , 
such that J{{x,y) = z x \/G{{x,y). The contours of the 
scalar stream function G{{x,y) = const correspond to 
streamlines of Jf (x, y), and we may chose Gf = for all 
points {x, y) along the inner and outer edges of the su- 
perconductor. Once a numerical result for <I>f is found, 
the result can be used to determine the effective areaSS, 
of the squid's central hole, A^s = <I>f/i?a, as was done 
for a circular ring in Ref. l3ll 

To prove that the effective area is also given by A^n = 
fnd/Id^ we consider the electromagnetic energy cross 
term = J{Bi ■ Bd//^o + A^oA^jf • jd)d.^r, where the 
integral extends over all space. Here, Bf{r) — Bg_{r) + 
Bsi{r) = V X A{{r), where j^{r) = V x Ba(r)/^o is 
the current density in the distant coil that produces a 
nearly uniform field Ba in the vicinity of the SQUID, 
jif = V X Bsi/no is the induced current density in the 
SQUID, and Bgi is the corresponding self-field under the 
conditions of flux focusing, i.e., when jif — through the 
junctions. Also, — W x A^ is the dipole-like held 
distribution generated by the circulating current Id with 
density jd in the SQUID; at large distances from the 
SQUIDSj^ Ad = A'o'Tid X T/47rr^. We evaluate £^fd in two 
ways, making use of the vector identities V • (A x B) — 
B • (V X A) - A • (V X B) and V • (7^) = 7V • j + V7 • j, 
and applying the divergence theorem, first with A = Ad, 
B = Bf, 7 = 7d, and j = jf, from which we obtain 



E{d = B^rnd, and then with A = Af, B = Bd, 7 = 7f, 
and j = jd, from which wc obtain ii^fd = ^fld with the 
help of Eq. ||2Jl. Since <I>f = i?a^cff, the effective area 
obeys A^s = md/h- 

The fourth term on the right-hand side of Eq. (|ll|l is 
due to a vortex at position — xxv + yyv in the body 
of the SQUID: 

*v(r-v) = / B^{x,y)dS + fioA ■ dl, (15) 

J s Jc 

where By{x,y) is the z component of the self-field gen- 
erated by the vortex's sheet-current density Jv{x,y) via 
the Biot-Savart law when no current flows through the 
junctions a and b. The desired fields are those that would 
appear in response to the vortex if the junctions a and 
b were open-circuited. Since V • Jv = 0, it is possible 
to express J^{x,y) in terms of a scalar stream function, 
as we did for Ji(x,y) and Jd{x,y). However, as shown 
below, it is possible to use energy arguments to express 
^v{x,y) in terms of the stream function Gd{x,y),^^ 

To obtain ^vi^) when a vortex is at the position 
r ~ XX + yy, imagine disconnecting the counterelectrode 
in Fig. 1 and attaching leads from a power supply to 
the contacts a and b. The power supply provides a con- 
stant current Id in the counterclockwise direction, and 
the sheet-current distribution through the body of the 
SQUID is given by Jd{x,y) = IdZ x \/Gd{x,y), as dis- 
cussed above. We also imagine attaching leads from a 
high-impedance voltmeter to the contacts a and b. If 
the vortex moves, the effective fiux <I>v changes with 
time, and the voltage read by the voltmeter will be^ 
Vab ~ d^^jdt. The power delivered by the power supply 
can be expressed in terms of the Lorentz force on the 
vortex, Jd X z^o = /d^oVGd; i.e., the rate at which work 
is done on the moving vortex is /d0oVGd ■ dr / dt. Equat- 
ing this to the power P = Idd^^/dt = /dV$v • dr/dt 
delivered by the power supply to maintain constant cur- 
rent, we obtain the equation \7^^{r) — 0oVGd(r). Thus 
<i>v(r) = (j)oGd{r)+const, where the constant can have 
one of two possible values depending upon whether the 
integration contour G is chosen to run inside or outside 
the vortex axis at r = xx+yy. Choosing G to run around 
the outer boundary of the SQUID, we obtain 

$v(r) = 0oGd(r). (16) 

Since Gd(To) ~ 1 for points r = ro on the outer edges 
of the SQUID and Gdivi) — for points r = on the 
inner edges (at the perimeter of the central hole or along 
the edges of the slit), we have ^v{ro) — <j>o and ^v{ri) — 
0. The derivation of Eq. (16 ) implicitly assumes that 
the vortex-core radius is much smaller than the linear 
dimensions of the SQUID. 

We now return to the problem of how to find the cur- 
rents I and Jd in the SQUID, as well as the critical 
current Ic- As discussed above, we have $/ = for a 
symmetric SQUID. For simplicity, we assume first that 
there are no vortices in the body of the SQUID, such that 
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FIG. 2: /c/27o vs <i>i/<f)o, calculated from Eqs. (O and itTHl . 

for TvLIo/cpo ~ 0, 1, 2, 3, 4, and 5 (bottom to top). 



<I)^ = and $ = + $d in Eq- ®7 where <I>d = Lid- 
From the sum and the difference of Ii and I2 we obtain 



I = 27o cos 



"0 



sin( 



Id 



In sin 



7r$f 7ri/d ' 

H — 1 cos < 

'0 00 



(17) 
(18) 



where (jj = (0i + (/'2)/2 is determined experimentaUy by 
how much current is appHed to the SQUID. When = 0, 
the current / is zero. As (j) increases, the magnitude of / 
increases and reaches its maximum value Ic at a value 
of (j) that must be determined by numerically solving 
Eqs. H17|l and H18|l . A simple solution is obtained for 
arbitrary $f only in the limit ttL/q/^o — s- 0, for which 
I = Ic = 2/0 1 cos(7r$t/0o)| and /d = at the critical cur- 
rent. For values of ttLIq/c/jq of order unity, as is the case 
for practical SQUIDs, one may obtain Ic for any value of 
$f by solving Eq. (|18|l self-consistently for /d for a series 
of values of and by substituting the results into Eq. H17|l 
to determine which value of maximizes /. Equations 
(|17|l and H18|) have been solved numerically by de Bruyn 
Ouboter and de Waele,^** (some of their results are also 
shown by Orlando and Dclin^^), who showed that at Ic 

icm = /c(*f + "0o) = /c(-$f), (19) 

/d(*f) = /d(*f + ri0o) = -/d(-$f), (20) 
/i($f) = /i($f + n0o) = /2(-$f), (21) 
hi^f) = /2($f + n0o) = /i(-$f), (22) 

where n is an integer. Hence all the physics is revealed 
by displaying Ic{^i) over the interval < $f < (j>o/2, as 
shown in Fig. 2. 

When a vortex is present, Eqs. H17|l and (|18|l still hold, 
except that <&£ in these equations is replaced by the sum 
$f + Thermally agitated motion of vortices in the 
body of the SQUID can produce flux noise via the term 
$v(^v) and the time dependence of the vortex position 



FIG. 3: Sketch of central portion of the long SQUID consid- 
ered in Sec. III. 



Tv. From Eq. we see that the sensitivity of Ic to 
vortex-position noise is proportional to the magnitude 
of V<I>v('") ~ 0oVGd = Jd X z(j)o/Id- Thus Ic is most 
sensitive to vortex-position noise when the vortices are 
close to the inner or outer edges of the SQUID, where 
the magnitude of Jd is largest. These equations provide 
more accurate results for the vortex-position sensitivity 
than the approximations given in Refs.|2i and<35i. 

So far, we have investigated how the general equations 
governing the behavior of a dc SQUID are altered when 
the contributions arising from line integrals of the current 
density are included. As we have shown in Ref. |^ these 
additional contributions are important when the Pearl 
length A is an appreciable fraction of the linear dimen- 
sions of the SQUID. We have found that the basic SQUID 
equations, Eqs. H17() and H18|l . are unaltered, except that 
the magnetic flux (sometimes called <I>exi~) generated in 
the squid's central hole by the externally applied fleld 
in the absence of a vortex is replaced by the effective 
flux $f , given in Eq. H14|l. Similarly, we have shown that 
the total inductance L of the SQUID has contributions 
both from the magnetic induction (geometric inductance) 
and the associated supercurrent (kinetic inductance). We 
also have shown in principle how to calculate the effect of 
the return flux from a vortex at position in the body 
of the SQUID, and we have found that the effective flux 
arising from the vortex is ^v{ry), given in Eqs. (|15|l and 
(|16|l . To demonstrate that all the above quantities can be 
calculated numerically for arbitrary values of A, we next 
examine the behavior of a model SQUID as decribed in 
Sec. III. 



III. LONG SQUID IN A PERPENDICULAR 
MAGNETIC FIELD 



Here we consider a long SQUID whose thickness d is 
less than the London penetration depth A and whose 
topology is like that of Fig. 1 but which is stretched 
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to a large length I in the y direction, as sketched in 
Fig. 3. SQUIDs of similar geom etry h ave been inves- 
tigated experimentally in Refs. l38(Tll87ll88ll39t We treat 
here the case for which the length I is much larger than 
the width 2w of the body of the SQUID, and we focus 
on the current and field distributions in and near the left 
{—w < X < —a) and right {a < x < w) arms and near 
the center of the SQUID, where to a good approximation 
the current density jy is uniform across the thickness 
and depends only upon x, and the magnetic induction 
B = V X A depends only upon x and z. In the equa- 
tions that follow, we deal with the sheet-current density, 
whose component in the y direction is Jy{x) = jy{x)d. 

The self- field magnetic induction generated by Jy {x) is 
B,j(x,z) = V X Aj(x,z), where Aj(x,z), the y compo- 
nent of the vector potential obtained from Ampere's law. 



Aj{x,z) = ^fjy{x')\u ^ 



The integration here and in the following equations is 
carried out only over the strips, and C is a constant with 
dimensions of length remaining to be determined. In the 
presence of a perpendicular applied field = zB^^ = 
V X Af, the total vector potential \s A — Aj + Ai, where 
= yB^x. 



A. Formal solutions 

We now use the approach of Ref.^^to calculate the in- 
plane magnetic-induction and sheet-current distributions 
appearing in Eqs. Hll(l - I|14|l in Sec. II. For all of these con- 
tributions we shall take into account the in-plane (z = 0) 
self-field contribution Aj{x) — Aj{x,Q) to the y compo- 
nent of the vector potential, where 



(24) 



We first examine the equal- current case and consider 
the contributions Bj{x), the z component of Bj{x), and 
Ji{x), the y component of Ji{x), due to equal cur- 
rents 1/2 in the left and right sides of the SQUID. 
Since Ji{~x) = Ji{x), the corresponding y component 
of the vector potential is also a symmetric function of 
x: Ai{~x) — Ai{x), where the subcripts / refer to 
the equal-current case. There are no flux quanta be- 
tween the strips (4>/ — 0), and the second term in 
the brackets on the right-hand side of Eq. Q vanishes; 
(0o/27i')V7 = 0. However, the constant C must be cho- 
sen such that Ji{x) = —Ai{x)/ (IqA in the superconduc- 
tor. Combining this equation with Eq. H24() . making 
use of the symmetry Jj{—x) = Ji(x), and noting that 
I = 2 Ji{x)dx, we obtain 



^ 1 ^ 



1 

2^ 



In- 



AS{x- 



Ji{x)dx' 
(25) 



for a < X < w. Here h can be chosen to be any conve- 
nient length, such as the length I or w, but not C. We 
now define the inverse integral kernel K^^{x,x') for the 
symmetric-current case via 

62 



K''yix,x") 



27r 



+ AS{x"-x') 



= Six- 

Applying this kernel to Eq. H25|) . we obtain 



Jiix) 



2tt 



■In 



Since 7 = 2/™ Jj{x)dx. 

.(^ 

such that 

/ 
2 



b 
C 

we find 



K''y{x,x')dx'. 



dx" 
x'). (26) 

(27) 



K''^{x,x')dxdx', 



(28) 



(23) Ji{x)^- K'y{x,x')dx'l K''^{x',x")dx'dx" 



(29) 



Gi{x) 



Ji{x')dx' , 



(30) 



X S -a, ^^i^-x) = -Gi{x). The 
z component of the magnetic induction 



For a < X < w, the stream function is 

2 

7 _ 

and for —w < x < —a, Gi{—x) 
corresponding 

Bi{x) can be obtained from the Biot-Savart law or, 
since Bi{x) — dAi{x)/dx, from Eq. (|^ . Note that 
Bi{-x) = -Bi{x). Although the kernel if^y (x, x') de- 
pends upon h, we find numerically that Ji{x) and Bj{x) 
are independent of b. 

We next examine the circulating- current case and con- 
sider the contributions Bd (x) , the z component of Bd {x) , 
and Jdix), the y component of Jd{x), due to a circulat- 
ing current Id; the current in the y direction on the right 
side of the SQUID is I2 — Id and that on the left side 
is /i = —Id- The vector potential is still given by Eq. 
(|^ . except that we add subscripts d. However, since 
Jd{—x) = —Jd{x), the vector potential is now an anti- 
symmetric function of x; Adi—x) — —Ad{x). The circu- 
lating current is generated by the fluxoid <I>d [see Eq. (|13|l ] 
associated with a nonvanishing gradient of the phase 7 
around the loop. The second term inside the bracket on 
the right-hand side of Eq. QJ, (0o/27r)V7, is -y<i>d/2; 
for a < a; < w and y$d/2Z for —w < x < —a. Thus, 
Jd{x) = -[Ad{x) - $d/2/]//ioA for a < a; < w. Com- 
bining this equation with that for Ad{x), noting that the 
inductance of the SQUID is L ^ 'I'dZ-^d, and making use 
of the symmetry Jd{~x) — —Jd{x), we obtain 



Lid 
21 



Mo 



27r 



X - 



A5{x-x') Jd{x')dx' (31) 



for a < X < w. We now define the inverse integral kernel 
K'^^{x^x') for the asymmetric-current case via 

1 



K'^^x.x"] 



2ti 



In 



A5{x" - x') 
= 5{x- 



dx" 



(32) 
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Applying this kernel to Eq. (|31|l and noting that Id 
Jd{x)dx, we obtain 



Jd{x) = aid I K''\x,x)dx' 



and 



where 



L = 2Q;/iQ/, 



a 1 / / / K''''{x,x')dxdx' 



(33) 



(34) 



(35) 



is a dimensionless function of a, 6, w, and A, which we 
calculate numerically in the next section. For a < x < w, 
the stream function is 



Gd{x) = ^ / Jd{x')dx', 



(36) 



and for —w < x < —a, Gd{—x) = Gd(x). 

In this formulation, as in Ref. L = + Ly, is 
the total inductance. The geometric inductance contri- 
bution = 2i?m//J, where -Em = I /™ Jd{x)Ad{x)dx is 
the stored magnetic energy, and the kinetic contribution 
Lk = 2£'k/-^d' where Ei^ — hqAI J^' J^{x)dx is the kinetic 
energy of the supercurrcnt, can be calculated using Eq. 
(123) from^-^l^ 



In 



X + x' 



Jd{x)Jd{x')dxdx' , 



ik = 



(37) 

^^°^^0|(.)d. = ^^, (38) 



(w - a) < Jd >2 



where the brackets (<>) denote averages over the film 
width. We can show that + L]^ ~ L with the help of 
Eqs. (1321) and 123) • 

The z component of the magnetic induction Bd{x) gen- 
erated by Jdix) can be obtained from the Biot-Savart 
law, or, since Bd{x) = dAd{x) /dx, from Eq. 124|) . Note 
that Bdi~x) = Bdix). 

When I ^ w, the magnetic moment in the z direction 
generated by the circulating current Id is (to lowest order 
in w/l) 



iTLd — 21 xJd{x)dx, 



(39) 



where the factor 2 accounts for the fac"b^ that the cur- 
rents along the y direction and those along the x direction 
at the ends (U-turn) give exactly the same contribution 
to nid, even in the limit Z — > cx). To next higher order in 
w/l one has to replace I in Eq. by I — {'w — a)q, where 
(ui — a)g/2 is the distance of the center of gravity of the 
x-component of the currents near each end from this end. 
For a single strip of width w — a, one has, e.g., (7=1/3 for 



the Bean critical state (with rectangular current stream 
lines) and q — 0.47 for ideal screening [K <^ w)^ 

We next examine flux focusing. As discussed in Sec. 
II, to calculate the effective area of the slot, we need to 
calculate the fields produced in response to a perpendic- 
ular applied magnetic induction _Ba, subject to the condi- 
tion that no current flows through either junction. Since 
this is equivalent to having both junctions open-circuited, 
the problem reduces to finding the fields produced in the 
vicinity of a pair of long superconducting strips connected 
by a superconducting link at only one end, i.e., when the 
slot of width 2a between the two strips is open at one 
end. However, the desired fields may be regarded as the 
superposition of the solutions of two separate problems 
when the slot has closed ends: (a) the fields generated 
in response to i?a, when <I> = (the zero-fluxoid case) 
and a clockwise screening current fiows around the slot 
[second term on the right-hand side of Eq. (|41|) below], 
and (b) the fields generated in the absence of i?a, when 
flux quanta in the amount of are in the slot and a 
counterclockwise screening current flows around the slot 
[first term on the right-hand side of Eq. 141|) ]. The de- 
sired flux-focusing solution is obtained by setting the net 
circulating current equal to zero. 

The equations describing the fields in the flux-focusing 
case are derived as follows. The z component of the 
net magnetic induction Bf{x) is the sum of i?a and the 
self- field Bs{{x) generated by Ji{x). The vector po- 
tential Ay{x) is the sum of Bi^x^ which describes the 
applied magnetic induction, and the self-field contribu- 
tion given by Eq. H24|) but with subscripts f. Since 
Jf{—x) = — Jf(a;), the vector potential is again an an- 
tisymmetric function of x; A{{--x) — —A{[x). The flux- 
oid $f [see Eq. (|14|l ] contributes a nonvanishing gra- 
dient of the phase 7 around the loop, such that the 
second term inside the brackets on the right-hand side 
of Eq. (QJ, ((/)o/27r)V7, is -y^(/2l ioT a < x < w 
and y^f/2l for —w < x < —a. Equation (Q yields 
Jf[x) = -[BaX + Af{x) ~ (f>f/2/]//xoA for a < a; < u;. 
Combining this equation with that for A{{x) [Eq. (|24() ]. 
making use of the symmetry J{{—x) = — Jf(x), and in- 
troducing the effective area via = i?a^cff [see Sec. II], 
we obtain 



Ba(Aeff - 2lx)/2l 



1^0 







X + x' 






'a 


V2tt 


X — x' 



+ A(5(a;-a;') Ji{x')dx' (40) 



for a < X < w. We again use the inverse integral kernel 
K^^{x,x') for the asymmetric-current case [Eq. (|32|) ] to 
obtain 

Mo J a " 2/ 

The effective area of the SQUID ^oft is found from the 
condition that the net current around the loop is zero 
[XT' Jii^)dx — 0], which yields 



Jf(x) 



{^-x')K^{x,x')dx'. (41) 



^eff — 2a/ / / xK^^{x, x )dxdx. 

J a J a 



(42) 
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For a < X < w, the stream function is 

Gf{x) = / Ji{x')dx', (43) 

J a 

and for —w < x < —a, G[{—x) — G{{x). The spatial 
distribution of the resulting z component of the in-plane 
magnetic induction is given by B{{x) = Ba + Bsf(x), 
where Bs{{x) can be obtained from the Biot-Savart law 
or by substituting J({x) into Eq. and making use of 
Bsi{x) = dA{{x)/dx. Note that B{{-x) = B{{x). The 
resulting magnetic moment rrif in the z direction can be 
calculated by replacing Jd by Jf in Eq. (jSHll- 

In the next section we also present numerical results 
for field and current distributions in the zero-fluxoid case, 
in which I = 0, = 0, and the effective flux is zero: 
$ = $d + $f = 0. Such a case could be achieved by short- 
circuiting the Josephson junctions in Fig. 1, cooling the 
device in zero field such that initially $ = 0, and then 
applying a small perpendicular magnetic induction i?a- 
A circulating current J{x), given by the second term on 
the right-hand side of Eq. 1)41(1 , would spontaneously arise 
in order to keep $ = 0, as in the Meissner state. 
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FIG. 4; Profiles of the sheet current Ji(x), Eq. I|49|l . stream 
function Gi{x), Eq. 15011 . and magnetic induction Bi{x) for 
the equal- current case (Ba = 0, Ii — I2 = //2 > 0). Shown 
are the examples a/w — 0.3 with A/w — (solid lines with 
dots), 0.03 (dot-dashed lines), 0.1 (dashed lines), 0.3 (dotted 
lines), and 1 (solid lines). Here B/no and J are in units h/w 
and G in units l2- 



where Sij = for i ^ j, Su — 1 , and 



B. Numerical solutions 

In the previous section we have presented formal so- 
lutions for the sheet-current density Jy{x) in Eqs. 
(|^ . and l|lTll . which are expressed as integrals involving 
the geometry-dependent inverse kernels K^'^{x,x') and 
K'^^lx, x'). As in Ref.i31. for thin rings, these integrals are 
evaluated on a grid (i = 1, 2, . . . , N) spanning only the 
strip (but avoiding the edges, where the integrand may 
have infinities), a < \xi\ < w, such that for any function 
f{x) one has f{x)dx = Wif{xi). Here Wi are the 

weights, approximately equal to the local spacing of the 
Xi\ the weights obey X^i^i Wi = w — a. We have chosen 
the grid such that the weights Wi are narrower and the 
grid points Xi more closely spaced near the edges a and w, 
where Jy{x) varies more rapidly. We have accomplished 
this by choosing some appropriate continuous function 
x(u) and an auxiliary discrete variable oc z — i, such 
that Wi = x'{ui){u2 — ui). We can choose x{u) such that 
its derivative x'{u) vanishes (or is reduced) at the strip 
edges to give a denser grid there. By choosing an appro- 
priate substitution function x{u) one can make the nu- 
merical error of this integration method arbitrarily small, 
decreasing rapidly with any desired negative power of the 
grid number N, e.g., N~'^ or N~^. 

For the equal-current case, Eq. (|25|l becomes 



16^ 



Mij)Ji{xj 



(44) 



27T \xf - Xj 



« 7^ J: 



2tt XjWj 



(45) 



The optimum choice of the diagonal term Q^f (i-e., with 
\xf — a;|| replaced by XiWi/ir for i = j, which reduces 
the numerical error from order N^^ to N^'^ or higher 
depending on the grid) is discussed in Eq. (3.12) of Ref. 
I41I for strips and in Eq. (18) of Ref. |^ for disks and 
rings. The superscript (sy) is a reminder that this is 
for a symmetric current distribution [Jj{~-x) — Ji{x)]. 
Defining K^J = {wjQ^j -\- Adij)~^ , such that 



N 
k=l 



and applying it to Eq. (|44|l . we obtain 



Jl{Xi) 



Since / — 2^^-^ WiJi{xi), we find 



N N 



such that 



1=1 j=l 



J N N N 



(46) 



(47) 



(48) 



(49) 



k=l 1=1 
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FIG. 5: Magnetic field fines in the equal- current case for 
a/w = 0.1 and K/w = 0, 0.01, 0.1, and 50 (or oo). 



It is remarkable that although the parameter b appears 
in Eq. the final result for Ji{xi) in Eq. ll^ does 

not depend upon b. The stream function Gi{x) can be 
evaluated aa^ 



i N N N 

j=i fc=i 



(50) 



1=1 m=l 



Shown in Fig. 4 are plots of Ji{x), Gi{x), and the cor- 
responding magnetic induction Bi{x) vs x for a/w = 0.3 
and various values of K/w = 0, 0.03, 0.1, 0.3, 1. The 
curves for A = exactly coincide with the analytic ex- 
pressions of Appendix A. The magnetic field lines for this 
case are depicted in Fig. 5. 

For the circulating- current case, Eq. (|31|) becomes 



21 



N 

Mo E(^j'3y' + A(5y) Jd(a;j) 



(51) 



where 



■In 



X 1 I X T 



27r \xi — 

1 A-KXi 

— In 

ZTT Wi 



(52) 



The superscript (as) is a reminder that this is for an 
asymmetric current distribution [Jd(— x) = — Jd(a;)]. 
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FIG. 6; Profiles of Jd{x), Eq. (|H|l, Gd{x), Eq. and mag- 
netic induction Bd{x) for the circulating-current case {Ba ~ 0, 
—Ii — I2 > 0). Shown are the examples a/w — 0.3 with 
A/w — (solid lines with dots), 0.03 (dot-dashed fines), 0.1 
(dashed lines), 0.3 (dotted lines), and 1 (solid lines). B/no 
and J are in units I2/U1 and G in units 12- 



Defining iff/ = {wjQf^ + A(5y)-i, such that 

N 



(53) 



fe=i 



applying it to Eq. (|51|l . and noting that 



Sill WiJd{xi), we obtain 



and 



where 



N N 

i/EE-^^^ 



(54) 



(55) 



(56) 



The stream function Gd (x) can be evaluated aa^ 



N 



Gd{x,) = aJ2J2^'3^7k■ 
J=l k=i 



(57) 



Shown in Fig. 6 are plots of Jd{x), Gd{x), and Bd{x) 
vs X for a/w — 0.3 and various values of A/w. Note that 
these curves look similar to those in Fig. 4, but they all 
have opposite parity, as can be seen from the different 
profiles B(x) near a; = 0. The magnetic field lines for 
this case are shown in Fig. 7. 

As discussed in Sec. II, when a vortex is present in 
the region a < \x\ < w, the sensitivity of the SQUID 's 



10 









FIG. 7: Magnetic field lines in the circulating -current case for 
a/w = 0.1 and A/w = 0, 0.01, 0.1, and 50 (or cxd). 



critical current Ic is proportional to the magnitude of 
d^v/dx — (f)i^dGd/dx = (f>oJd{x)/ Id- From Fig. 6 we see 
that when A <C ic, this sensitivity is greatly enhanced 
when the vortex is close to the edges a and w but that 
when A > w, the sensitivity is nearly independent of 
position. 

Shown as the solid curves in Fig. 8(a) are plots of the 
inductance L vs a/w for various values of A/w = 0, 0.03, 
0.1, 0.3, and 1. The solid curves in Fig. 9(a) show the 
same L vs A/w (range 0.0045 to 2.2) for several values of 
a/w = 0.01, 0.1, 0.4, 0.8, 0.95, and 0.99. 

The geometric and kinetic contributions and 
can be calculated separately from Eqs. (|37|) and H38|) 



ini = ^^^y^^y^^WiWjQ'^Jd{xi)Jdixj) , (58) 

(59) 



N N 



1=1 j=l 
N N 



ik = '^^J^^^W^Jj{x,) , 

d i=l j=l 



using Eqs. H54|) and H56|l . We can show that Li„ + Lk = L 
using the property of inverse matrices that M ■ M^^ = 
■ M = J, where I is the identity matrix. 
Shown as solid curves in Fig. 8(b) are Lm and Lk vs 
a/w. For A = 0, when ik = 0, L = exactly co- 
incides with Eq. (AlO), which may be approximated by 
Eq. (All) for a/w < 0.7 [open circles in Fig. 8(a)] and 
by Eq. (A12) for a/w > 0.7 [open squares in Fig. 8(a)]. 
The dotted curves in Fig. 8(b) for L,„ and Lk are those 
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FIG. 8: (a) Solid curves show the inductance L — Lm + ik 
[Eq. (Hg] vs a/w calculated for A/w = 0, 0.03, 0.1, 0.3, and 
1. See the text for descriptions of analytic approximations 
shown by the open circles, open squares, dots, and dashes, 
(b) Solid curves show the geometric inductance Lm [Eq. 1581 1 
vs a/w for A/w = 0, 0.1, and 1 and the kinetic inductance 
Lk [Eq. lEU] vs a/w for A/w = 0.03, 0.1, 0.3, and 1. See the 
text for descriptions of analytic approximations shown by the 
dotted and dashed curves. 



of Eqs. (C3) and (C6) in the limit A/w oo, when the 
circulating current density is uniform. The dotted curves 
in Fig. 8(a) are obtained from L — + Lk using the 
approximations of Eqs. (C3) and (C6); they are an excel- 
lent approximation to L for A/w > 0.03 except for small 
values of a/w. Improved agreement for small values of 
A/w and a/w is shown by the dashed curves in Fig. 8(a), 
which show the approximation of Eq. (B7) for L, and in 
Fig. 8(b), which show the approximation of Eq. (B12) for 
Lk. 

The solid curves in Fig. 9(b) show i,,! and ik vs A/w. 
The geometric inductance depends upon A but only 
weakly, varying slowly between its A = asymptote 
[Eq. (AlO), horizontal dot-dashed line] and its A = oo 
asymptote [Eq. (C3), horizontal dotted line]. For larger 
values of of a/w, Lm is nearly independent of A. On 
the other hand, the kinetic inductance L^, is approxi- 



11 




- 10" 

o 
zi. 



- (b) ' x^L, 

a / w =/0.99 


1 






- / ^ 


/0.95 






/ / /0.8 








m ; 


■ / ^^^^L—^^^- 








L 

m ■ 


^ 1 1 1 1 1 1 1 1 


1 


1 1 1 


, , , , 1 



10"' 



10" 
A/ W 



10" 



FIG. 9: (a) Solid curves show the inductance L — Lm + ik 
[Eq. JSSJ] vs A/w for a/w = 0.01, 0.1, 0.4, 0.8, 0.95, and 
0.99. See the text for descriptions of analytic approximations 
shown by the dotted and dashed curves, (b) Solid curves show 
both the geometric inductance Lm [Eq. ((HHl] a-nd the kinetic 
inductance Lk [Eq. lIH^ I vs A/w for a/w = 0.01, 0.1, 0.4, 
0.8, 0.95, and 0.99. See the text for descriptions of analytic 
approximations shown by the dotted, dashed, and dot-dashed 
curves. 
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FIG. 10: The magnetic moment for the circulating current 
case [Ba ~ 0, — Ji ~ I2 ~ Id > 0) plotted versus a/w for 
various values of A/w in units 2wlld- These curves coincide 
with those in Fig. 13 below. 
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mately proportional to A/w. The straight dotted hnes 
in Fig. 9(b), calculated from the large-A approximation 
given in Eq. (C6), are a good approximation to L^^ ex- 
cept for small values of A/w and a/w. The dotted curves 
in Fig. 9(a) are obtained from L = Lm + using the 
approximations of Eqs. (C3) and (C6). Improved agree- 
ment for small values ol A/w and a/w is shown by the 
dashed curves in Fig. 9(a), which show the approxima- 
tion of Eq. (B7) for L, and in Fig. 9(b), which show the 
approximation of Eq. (B12) for Lk- 

In Eq. (3) of Ref. |3 Yoshida et al. derived an ap- 
proximate expression for the kinetic inductance when 
A/w ^ 1. We have found that their expression for L-^ 
is not an accurate approximation to our exact numer- 
ical results. To eliminate the logarithmic divergences 
due to the inverse square-root dependence of the cur- 
rent density near the edges, Yoshida et al. followed an 



FIG. 11: Profiles Jf{x), Eq. JHSJ, Gf{x), Eq. igij, and mag- 
netic induction Bf{x) for the flux-focusing case [Ba > 0, 
Ii = I2 ~ 0). Shown are the examples a/w = 0.3 with 
A/w — (solid lines with dots), 0.03 (dot-dashed lines), 0.1 
(dashed lines), 0.3 (dotted lines), and 1 (solid lines). B and 
fj,oJ are in units Ba, and G in units wBa/jJ^o- 



approach used by Meservey and Tedrow^a^ and chose a 
cutoff length of the order of d, the film thickness. When 
d < A, however, this approach cannot be correct, be- 
cause the equations describing the fields and currents in 
superconducting strips contain only the two-dimensional 
screening length A — X^/d. The cutoff length therefore 
must instead be chosen to be of the order of A, as we 
have done in Appendix B. 

The magnetic moment in the z direction generated by 
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the circulating current is, from Eq. (|39|l . 

N 

md = 2 ^ WiXiJdixi). (60) 

i=l 

As shown in Fig. 10, this magnetic moment vanishes very 
slowly when the gap width and A go to zero, a/w — s- 
and A/w — > 0. This can be explained by the fact that for 
A = and a < x w one has Jd{x) oc 1/x, Eq. (A5). 
The contribution of these small x to nid, Eq. H60() . stays 
finite due to the factor x, but the total current Id to 
which nid is normalized, diverges when a/w —^ 0, thus 
suppressing the plotted ratio irid/Id- Interestingly, the 
curves in Fig. 10 coincide with those in Fig. 13; see below. 
Expressions for in the limits A/w — > and A/w oo 
are given in Eqs. (A13) and (C7). 

For the flux-focusing case, Eq. 140|l becomes 

N 

B,{A,s/2l - X,) = Y^^WjQTj + A5,j)Ji{xj). (61) 
i=i 

Applying Eq. lf5S|) . we obtain 
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FIG. 13: The effective area Acff , Eq. I63II . plotted versus the 
gap half width a/w for several values of A/w — 0, 0.001, 0.003, 
0.01, 0.03, 0.1, 0.3, and 10. The lower-right curves show the 
same data shifted and stretched along a/w. The dots show 
the exact result (A19) in the limit A/w 0. For A/a > 10 
one has Ac:ft/2wl ~ (1 + a/w)/2, Appendix C. These curves 
coincide with Fig. 10 since A^s — nid/ Id- 




A/w 

FIG. 14: The effective area A^s, Eq. (1631 . plotted versus A/w 
(range 7 • 10~^ to 14) for several values of a/w = 0, 0.001, 
0.003, 0.01, 0.03, 0.1, and 0.2. Same data as in Fig. 13. 

where, since X^i^i WiJf{xi) = 0, the effective area is 

N N 

Aeff = 2aZ ^ ^ w.K^Xj . (63) 
The stream function Gi{x) can be evaluated aa^ 

i 

Gf(x,) = ^u;,Jf(x,). (64) 
i=i 

Shown in Fig. 11 are plots of the flux- focusing Ji{x), 
Gf{x), and B{{x) vs x for a/w — 0.3 and A/w = 0, 0.03, 
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FIG. 15: The minimum of the magnetic induction in the flux- 
focusing case, B{{0) = B{{x = 0), referred to the applied field 
Ba and plotted versus the half gap width a/w. Top: The 
ratio B{{x = 0)/Ba, tending to unity for a/w — > 1 and for 
A/w ^ 1, and diverging for a/w — > when A = 0. Bottom: 
The same ratio multiplied by a/w to avoid this divergence 
and fit all data into one plot. The lower right plot depicts 
the small-gap data two times enlarged along the ordinate, 
and shifted and five times stretched along the abscissa. The 
circles show the approximation Bf{0)/Ba ~ (w/a)/ ln(4w/a) 
good for a/w < 0.34^ 



0.1, 0.3, and 1. The first term in Eq. (|62|l equals the 
circulating-current sheet-current density, Eq. (|54|l . with 
appropriate weight factor such that the total circulating 
current vanishes, Ii — I2 — 0. The corresponding mag- 
netic field lines are depicted in Fig. 12. Shown in Figs. 13 
and 14 are plots of the effective area Aefi{a/w, A/w) ver- 
sus a/w and A/w, respectively, in units of the maximum 
possible area 2wl. In the limit A/w — *■ 0, A^s is given 
by Eq. (A19), and when A/w 00, A^ft = l{w + a). 
Note in Fig. 14 that Acs increases with increasing A, 
particularly for small gap widths 2a. Flux focusing is re- 
flected by the fact that for small a/w ^ the effective 
area A^s of the gap tends to a constant, except in the 
limit A — > 0, where it vanishes very slowly, Aef[/2wl ~ 
(7r/2)/ln(4-u;/a) [Eq. (A19)]. When a/w 0, the en- 
hancement factor Aeff/2al — s- 00 and thus diverges even 
for A = 0. In the limit a/w ^ 1, Acs{0,A/w) tends to 




A / w 

FIG. 16: The minimum field Bt{0) = B{(x = 0) for the flux- 
focusing case as in Fig. 15 but plotted versus A/w (range 
7 ■ 10"® to 14) for several values of a/w = 0, 0.001, 0.003, 
0.01, 0.03, 0.1, and 0.2. 
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FIG. 17: The magnetic moment m/ for the flux-focusing 
case plotted versus a/w for various values of A/w in units 
w'^lBa/fJ.o- 



a universal function [see Fig. 14]. Interestingly, Figs. 13 
and 10 show identical curves; this is because the identity 
Acff — rud/Id holds for all values of a/w and A/w, as 
proved in general in Sec. II. 

Figures 15 and 16 show the minimum of the magnetic 
induction in the flux-focusing case, Bf{0) — B{{x = 0) 
[see Fig. 11], plotted versus a/w and A/w, respectively. 
The ratio B[{0)/Ba > 1 tends to unity for a/w 1 
and for A/w 3> 1, and it diverges for a/w —^ when 
A = 0. The curve for A = exactly coincides with 
the analytic expression Bf{0)/Ba = wE{k')/aK{k') 
obtained from Eq. (A17). For a/w ^ 1 this yields 
Bf{0)/Ba ~ {w / a) / \n{Aw / a) , which is a good approx- 
imation for < a/w < 0.3.-'^^ The magnetic moment nif 
for the flux-focusing case, calculated from Eq. H60|l but 
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FIG. 18: Profiles of tlie sfieet-current density J(x) [second 
term in Eq. 1621 1 and magnetic induction B{pc) generated 
wfien a perpendicular magnetic induction Bz. is applied in 
the zero-fluxoid case when $ = 0, 7i = —I2 > 0, and / = 0. 
Shown are the examples a/w = 0.3 with A/w = (solid lines 
with dots), 0.03 (dot-dashed lines), 0.1 (dashed lines), 0.3 
(dotted lines), and 1 (solid lines). B and fioJ are in units Ba- 



with Jd{xi) replaced by Ji{xi), is shown in Fig. 17. Ex- 
pressions for TOf in the limits A/w — *■ and A/w — + 00 
are given in Eqs. (A20) and (C9) 

Shown in Fig. 18 are profiles for the zero-fluxoid case 
with plots of J{x) and the corresponding B{x) gener- 
ated by an applied magnetic induction Ba > when 
the junctions are short-circuited such that $ = and 
1 1 — —I2 > 0; for comparison see analogous profiles in 
Sec. 2.5 of Ref. for two parallel strips and in Sec. 
IV of Ref. m and Sec. 4 of Ref. iH for rings. That the 
current density J{x) in the zero-fluxoid case is given by 
the second term on the right-hand sides of Eqs. H41|l and 
(|H2J), can be seen by setting $f = B^Acs = in Eqs. 
601, (EJ), (EH), and (EH. Depicted in Fig. 18 are the 
examples a/w = 0.3 with A/w = 0, 0.03, 0.1, 0.3, and 
1. Figure 19 shows the magnetic field lines for this case 
and Fig. 20 the magnetic moment m. Expressions for J, 
B, and m for the zero-fluxoid case in the limits A/w — s- 
and A/w ^ 00 are given in Appendixes A and C. 



IV. SUMMARY 

In Sec. II of this paper we have presented general 
equations governing the static behavior of a thin-film dc 
SQUID for all values of the Pearl length A = /d, where 
the London penetration depth A is larger than d, the 
film thickness. The SQUID's critical current Ic depends 
upon the effective flux <&, which is the sum of the mag- 
netic flux up through a contour surrounding the central 
hole and a term proportional to the line integral of the 
current density around this contour. For a symmetric 
SQUID there are three important contributions to <I>: a 
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a/w 

FIG. 20: The magnetic moment m for the zero-fluxoid case 
^ — 0, Ba > 0, and 1 = plotted versus a/w for various 
values of A/w in units w^lBa/ fJ-o- At a = A = one has 

— m = TTW^lBa/ fJ-O- 



circulating-current term <I>d, a vortex-field term <I>v, and a 
fiux-focusing term <I>f , all of which depend upon A. Since 
A is a function of temperature, an important consequence 
is that all of the contributions to are temperature- 
dependent. 

The circulating-current term $d can be expressed in 
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terms of the SQUID inductance L and the circulating 
current via $d — LI^. The SQUID inductance has 
two contributions, L = + Lk, where the first term 
is the geometric inductance (associated with the energy 
stored in the magnetic field) and the second is the kinetic 
inductance (associated with the kinetic energy of the cir- 
culating supercurrent). Both contributions are functions 
of A, since they both depend on the spatial distribution 
of the current density. However, depends only weakly 
upon A, because for the same circulating current 1^, the 
energy stored in the magnetic field does not vary greatly 
as A ranges from zero to infinity. On the other hand, 
because the kinetic energy density is proportional to A, 
Lk is also nearly proportional to A, with deviations from 
linearity occurring only for small values of A/w. 

The vortex- field term can be written as $v = ^oGd, 
where Gd is a dimensionless stream function describing 
the circulating sheet-current density Jd- Roughly speak- 
ing, when A is small, Ic is most strongly dependent upon 
the vortex position when the vortex is close to the edges 
of the film, but when A is large, Ic is equally sensitive 
to the vortex position wherever the vortex is. Recent 
experiments^ have used the relationship $v = 0oGd to 
determine the vortex- free sheet-current density Jd(x,y) 
from vortex images obtained via low-temperature scan- 
ning electron microscopy.^'2> The experimental data ob- 
tained in magnetic fields up to 40 /iT are in excellent 
agreement with numerical calculations of Jd{x,y), con- 
firming the validity of the above relationship, even in 
the presence of many (up to 200) vortices in the SQUID 
washer. 

The flux-focusing term can be expressed as $f = 
_Ba^cff, where i?a is the applied magnetic induction and 
Acs is the effective area of the central hole of the SQUID. 
Although Acs is primarily determined by the dimensions 
of the SQUID, it also depends upon the value of A. 

To illustrate the A dependence of the above quantities, 
in Sec. Ill of this paper we analyzed in detail the be- 
havior of a long SQUID whose central region resembles a 
coplanar stripline. We numerically calculated the profiles 
of the sheet-current density, stream function, and mag- 
netic induction in the equal-current, circulating-current, 
flux-focusing, and zero-fluxoid cases for various represen- 
tative values of A. We presented plots of the inductances 
i, im, and Lk, the effective area ^eff, and the magnetic 
moments for these cases. Useful analytic approximations 
are provided for the A/w ^ limit in Appendix A, for 
small A/w and a/w in Appendix B, and for the A/w — > oo 
limit in Appendix C. 

We are in the process of applying the above theory to 
square and circular SQUIDs, using the numerical method 
of Ref.EO- 
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APPENDIX A: THE LIMIT A/w = 

In the ideal-screening limit A/w = 0, the y component 
of the sheet-current density in the strips {a < \x\ < w) 
for the equal- current case is^ 



Jiix) 



7r[(.T2-a2)(ti;2-x2)]i/2' 



(Al) 



and the z component of the magnetic induction in the 
plane z = of the strips is 



Bi{x) 



tiol 



2n [{x^ - a^){x^ - w^)]^/ 
0, a < < w, 
HqI X 



|a;| > w, (A2) 
(A3) 



2tt [(a2_:^2)(^2_3.2)]i/2' 1^1 (A4) 

and the constant C in Eqs. ^7^, and ijlHll is 

C = Vu'2 -a2/2. 

For the circulating- current case, the y component of 
the sheet-current density in the strips (a < |a;| < w) in 
the hmit A/w ^ ig& 



Jd{x) = 



2Ba x 



fio \x\ [(a;2-a2)(w2-a;2)]i/2' 



(A5) 



and the z component of the magnetic induction in the 
plane z = of the strips is 



w 



Bdix) = -Boj-^ 1^1 (^^) 

= 0, a < |a;| < w, (A7) 



Bo 



[(a2 -2;2)(w2 _2;2)]1/: 



— , \x\ < a, (A8) 



where the parameter Bq, the magnetic flux (f>d in the z 
direction in the slot, the circulating current Id, and the 
geometric inductance Lm are related by 



and 



*d = L^Id = 2BolwK{k) 



- tiolK{k)/K{k'), 



(A9) 



(AlO) 



where K{k) is the complete elliptic integral of the first 
kind of modulus k = a/w and complementary modulus 
k' = Vl — fc^. The geometric inductance is well approx- 
imated for small a/w by 



im = (7r^oV2)/ln(4w/a), 



(All) 
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neglecting corrections proportional to a^/uP', and for 
small {w — a)/w by 

Lrn = {^iol/7^) ln[16/(l - a^w^ (A12) 

neglecting corrections proportional to 1 — a^/w'^. In the 
limit that A = 0, the kinetic inductance vanishes (L^ — 
0), and the inductance in Eq. (|A10p becomes the total 
inductance: L = im. The magnetic moment [see Eq. 
(EnH can be obtained from Eqs. (13)-(16) of Ref.lil 

TOd = [nlw/K{k')]Id. (A13) 

For the flux-focusing case, the y component of the 
sheet-current density in the strips (a < |a;| < w) in the 
limit A/w = is'*'^ 



Jf(x) = 



2Ba X E{k')w'^ - 2K{k')x 



2 -a2)(y,2 _2;2)li/2' ^^^"^^ 



^ioK{k') \x\ [(a;2 -a2)(u;2 -2;2)] 



and the z component of the magnetic induction in the 
plane 2 = of the strips is 

. , E{k')w^ - 2K{k')x^ 

^^•^^ K{k')[{x^ -a^){x^ -w^)]^/^' 

\x\ > w, (A15) 

= 0, a < |a;| < w, (A16) 

E{k')w^ - 2K(k')x^ 
K{k') [(a2 -a;2)(u;2 _a;2)]i/2' 

< a, (A17) 

where E{k') is the complete elliptic integral of the sec- 
ond kind of complementary modulus k' = y/l — k^ and 
modulus k ~ a/w. The magnetic flux in the z direction 
in the slot is 

i^>f = TrBJw/K{k'), (A18) 
and the effective area Acff of the slot is 

Aeff = $f/Ba = TTlw/K{k'). (A19) 

Note that Acs = nid/Id- The magnetic moment gener- 
ated by Ji{x) is 

mf = -nl[w^ - 2u? E{k') / K{k')]B^/ ii^. (A20) 

For the zero-fluxoid case, the y component of the sheet- 
current density in the strips can be obtained from Sec. 
2.5 of Ref. 47: 

2B, X x^-[l-E{k)/K{k)]w^ 
^ ' Mo \x\ [(x2-a2)(^2_a.2)]i/2- ^ ) 

The corresponding z component of the magnetic induc- 
tion isii 

x^-\l-E(k)/K(k)]w^ , , , 



[(a;2-a2)(^2_^2)]i/2 
= 0, a < |x| < w, 



(A23) 



_ [l^E(k)/K{k)]w^-x^ 

^ ^ [(a2-a;2)(^2_^2)]i/2'l^l <«• 

The magnetic moment generated by J{x) is 

m = -■Kl[2w'^E{k)/K{k) -w"^ + a^]Bjfio. (A25) 



APPENDIX B: BEHAVIOR FOR SMALL A AND 
SMALL a 

In this section we present some expressions for L, L]^, 
and Lin that follow from approximating the circulating- 
current distribution for small values of A and a. 

When the slot is very narrow [a/w <C 1), we approxi- 
mate the sheet-current density in the region a < x < w 
generated by the fluxoid $d via 

w 

Jdix) = /q , (Bl) 
V(a;'-a^+'52)(ii'2_^2 + j2)' 

where /q = 2<I>d/7r/xo? and (5 is a quantity of order A = 
\^ /d determined as follows. When I 00 and then a 
and w 00, an exact calculation yields for a; > 



Jy{x) 



2$d f°° e 



(B2) 



We find J^Jyix)dx = (2$d/7rAioO M7V2A) when b > 
A, where 7 = e*^ = 1.781..., and C = 0.577... is Eu- 
ler's constant. From Eq. I|B1|) we find fl^ Jd{x)dx = 
(2<i>d/7rMoO ln(26/5) when a = and i5 < 5 < 
Comparing these two integrals we obtain S — (4/7) A — 
2.246A. Integrating Eq. ljBl|) from a to w to obtain 1^, 
we find 

h = Io^^==[F{Xa,q)-F{\^,q)], (B3) 
where k) is the elliptic integral of the first kind and 



Aq = arcsin 
A7/1 = ai 



<52 



2<52 ' 



y/up 



2,52 



252 



(B4) 
(B5) 

(B6) 



Expanding Eq. (|B3|) for a <C w and (5 ^ w, using 
/q = 2<I>d/7r/xoZ, and neglecting terms of order a^ /w^ and 
juP' , we obtain 

i = $d//d = (^/^oV2)/{ln[4u;/(a -I- <5)] - bjw}, (B7) 

where 5 = 2.246A. Note that Eq. (B7) reduces to Eq. 
(All) when A = 0. 
From Eq. ((BlJ we obtain the approximation 



j\{x)dx — /q 



w 



(52) 



(/a + /^), (B8) 



where 



j a 



W 



: tan 



_i (w — a)\/<52 — a2 
a(u' — a) + (52 



a < 6, 
w 



(B9) 



- ,52 
a > (5, 

Vw^ + (52 



: tanh' 



(it; — a)y/a'^ — 6^ 
a{w — a) -I- (52 ' 



(BIO) 

1 {w-a)Vw^+S'^ 

tanh rT^2~- (^11 

w(w — a) -I- 
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Using Eqs. (jBSp and (|B8p . we obtain from Eq. 

.A (W^+S^) ifa + U) 



w {w^ -a^ + 2(52) [F{X,,q) - F(A^, g)]2 ' 

where (5 = 2.246A. Although our intention in using the 
ansatz of Eq. I|B1I) initiaUy was to obtain an improved 
approximation to Lk for small values of A and a, we 
see from Figs. 8(b) and 9(b) that Eq. (|B12|) provides a 
reasonably good approximation for all values of A and a. 



APPENDIX C: THE LIMIT A/w ^ oo 

In the weak-screening limit A/w —>■ oo, K^^(x,x') = 
A~^d{x — x'), the y component of the sheet-current den- 
sity in the strips (a < < w) in the equal-current case 
is uniform, J/ = I /2{'w~a), the z component of the mag- 
netic induction in the plane of the strips obtained from 
the Biot-Savart law is 



Biix) 



IJ-oI 



■In 



(x ~ w){x -\- a) 



4:Tt(w — a) (x + w){x ~ a) 



(CI) 



and the constant C in Eqs. (EH), (ESJ, and ^ is 

C — w exp[— 7rA/(?ii — a)]. 

For the circulating- current case in the limit A/w — > 
oo, K''''{x,x') = A~i(5(a; - x'), a = A/{w - a), the y 
component of the sheet-current density in the strips is 
again uniform, Jd = Id/ {w ~ a) for a < x < w, and the z 
component of the magnetic induction in the plane of the 
strips obtained from the Biot-Savart law is 



Bd{x) 



tiold 



■In 



w 



2tt{w — a) 

The geometric inductance is, from Eq. H37|) 



(C2) 



^lQl 



tt{w — a)2 
—2aw In 



w In 



\w — a'^ J 



In 



4a2 



w^ 



(C3) 



which is independent of A. Equation (C3) is well approx- 
imated for small a/w by 



Ln, = {nol/TT){l + 2a/w) In 4, 

neglecting corrections proportional to a^/w"^, 
small {w — a)/w by 



In 



1 — a/w 



- - -(I -a/w) 
2 2^ ' ' 



(C4) 
and for 

, (C5) 



neglecting corrections proportional to (1 — a/w)'^. From 
Eq. H38(l we obtain the kinetic inductance 



Ly_ ~ 2/io^A/ (w ~ a). 



(C6) 



When A^ w, the total inductance L is dominated by the 
kinetic inductance (Lk 3> im), such that i w ik- Since 
Jd is uniform, the magnetic moment is easily found from 
Eq. ^ to be 



md = l{w -\- a)/d. 



(C7) 



oo, 

the y component of the sheet- 
= B^{w -\- a - 2x)/2plqA for 



For the flux-focusing case in the limit A/w 
K'^^x.x') = A-^5{x-x 
current density is Jf{x) 
a < X < w, the effective area is Acb = l{w + a) = rrid/Ii, 
and the z component of the magnetic induction in the 
plane of the strips is B{{x) = Bi^ -\- Bsi{x), where from 
the Biot-Savart law 



Bsi{x) 



B, 



47rA 
+ 2xln 





X^ — W^ 


{w + a) In 


? 9 

X — a 



{x + w){x — a) 



■i{w - a) . (C8) 



{x — w){x -\- a) 
The magnetic moment generated by J{{x) in this limit is 



mt = ~[l{w-^af/6A]{Bjfio). 



(C9) 



For the zero-fluxoid case in the limit A/w — > oo, the 
applied field is only weakly screened, and the z compo- 
nent of the magnetic flux density is nearly equal to the 
applied magnetic induction, B{x) « Bg_. The y compo- 
nent of the vector potential is approximately given by 
A{x) = Bi^x, and the y component of the induced sheet- 
current density, obtained from Eq. with 7 = 0, is 
J{x) — —{Bs,/ noA)x. To the next order of approxima- 
tion, B{x) = Bs. + Bs{x), where the self- field Bg is found 
from the Biot-Savart law 



Bsix) = 



B, 
2nA 



:ln 



{x + w){x — a) 



-2{w-a) . (CIO) 



(a; ~ w){x -\- a) 
The magnetic moment generated by J{x) in this limit is 
m = -[2l{w^-a^)/3A]iBjfio)- (CU) 
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